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I. INTRODUCTION 


The electromagnetic radiation from lightning at microwave frequencies and below is generally referred to as 
the “radio frequency” or “RF” portion of the spectrum. Radiation from lightning in this portion of the spectrum 
is important both for scientific investigations of lightning and for engineering assessments of the interference envi- 
ronment during thunderstorms. Measurements have been reported from frequencies below a kilohertz to frequen- 
cies above a gigahertz. 

The word “spectrum” is generally used in the literature on lightning in this frequency range to mean the 
magnitude of the Fourier transform of the electric field E(t) radiated during the discharge. In applications to light- 
ning, two methods have traditionally been employed to measure this spectrum. In one, the spectrum is obtained 
from the electric field waveform itself by Fourier transforming. The electric field waveform is recorded first using 
instrumentation such as a fast field change system and wide bandwidth recorder, and the spectrum is obtained 
numerically from this record by means of a Fourier transform. This technique has the appeal of being straight for- 
ward, but it requires wide bandwidth recording devices with large dynamic range because the power at high fre- 
quencies tends to decrease rapidly with frequency. The second technique is to measure the energy radiated at a 
particular frequency directly using a filter and detector system tuned to the frequency of interest. Standard radio 
receivers suitable for this purpose are available in this frequency range. The major difficulty with measurements 
of this type has been in identifying the element of the lightning flash (e.g., leader, return stroke, etc.) which is the 
source of the radiation. Each of these techniques will be described in this report and a summary given of the data 
obtained with each. 


II. SPECTRA OBTAINED BY FOURIER TRANSFORM 

A lightning flash is not a single event but, rather, is a sequence of many discrete events. Some, such as return 
strokes and the leader steps preceeding a first return stroke, have received much attention and are reasonably well 
known, whereas others such as the many different events which take place in the cloud are only recently beginning 
to be studied. The electric fields radiated from these events tend to have characteristic (time domain) shapes which 
permit the events to be identified and spectra to be obtained separately for the various events. 

For example, at the top in Figure 1 is shown an electric field waveform E(t) recorded by the author in Florida 
during the Thunderstorm Research International Project (TRIP; Pierce, 1976). This waveform has a shape charac- 
teristic of those observed during first return strokes (Uman and Krider, 1982; Weidman and Krider, 1978). It be- 
gins abruptly with a rapid rise to peak and then decays irregularly toward zero. Frequently, the abrupt beginning is 
preceeded by a string of small pulses associated with the stepped leader (Krider et al. , 1977; Weidman and Krider, 
1980). The last such step is visible in this record just before the beginning of the return stroke. The graph at the 
bottom of Figure 1 is the magnitude of the Fourier transform of the waveform, E(t), shown at the top. The Fourier 
transform was obtained numerically after the waveform at the top was digitized. 

Figure 2 shows the results of averaging waveforms from 20 return strokes. The data was recorded in Florida 
during a short period of strong activity associated with a nearby thunderstorm on July 20, 1976. At the top is the 
average of the time domain electric field waveforms, E(t), and at the bottom is the average of the magnitude of the 
Fourier transform (the magnitudes were averaged). This data was collected by the author during TRIP-76. 

One of the earliest applications of the Fourier transform approach was made by Watt and Maxwell (1957). In 
the course of describing measurements of low frequency attenuation due to the earth, they made Fourier transforms 
of electric field waveforms recorded by Norinder (1954) and Florman (1955) and from them formed a composite 
spectrum; however few details were given regarding the waveforms. Later, Taylor (1963) reported a carefully 
documented measurement of the spectrum in which the lightning discharges were located and propagation losses 
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Figure 2. Average electric field waveform (top) and average magnitude of the Fourier 
transform of the individual waveforms (bottom). Data recorded from first return 
strokes in Florida during TRIP-76. 







taken into account. Taylor recorded return stroke waveforms and made Fourier transforms over a frequency range 
from 1-100 kHz. The spectrum he obtained is shown in Figure 3 (open circles) normalized to a range of 50 km. 
The normalization has been done assuming a (distance)' 1 dependence for the amplitude. The spectrum obtained by 
Taylor peaks near 5 kHz and then decays roughly as (frequency)' 1 to 100 kHz, the limit of the measurements. In 
obtaining this spectrum, Taylor used the recently verified theory for the effects of attenuation due to propagation 
over the earth to correct for the loss of signal at low frequencies. He had both this theory and tools for locating the 
lightning discharge at his disposal. 

The Fourier transform approach was largely ignored until Serhan et. al. , (1980) again reported spectra of radi- 
ation from return strokes. Using modem techniques for locating nearby return strokes (Krider et. al., 1976; Krider 
et. al., 1980), Serhan et. al. were able to separate first and subsequent return stroke waveforms and compute 
spectra for each. Their data for first return strokes at 50 km are shown in Figure 3. These data fit well with Taylor’s 
(1963) measurements and extend the (frequency)' 1 trend to several hundred kHz. The spectra obtained by Serhan 
et. al. (1980) for subsequent return strokes has substantially the same shape but is somewhat lower in amplitude. 

Weidman et. al., (1981) using substantially similar techniques also reported spectra of first return strokes 
(squares in Figure 3). Their data, collected from lightning over the ocean to minimize effects of propagation, ex- 
tend the spectrum to 1 MHz and continue to show a decrease proportional to (frequency)" 1 . 

Weidman et. al., (1981) also reported first return stroke spectra at higher frequencies which they obtained by 
recording the derivative, dE/dt, of the electric field waveform rather than E(t) itself. Recording the derivative im- 
proves the sensitivity of the measurement to higher frequencies (because the spectrum of the derivative is the spec- 
trum of the waveform multiplied by frequency). The data obtained by Weidman et. al., (1981) in this manner are 
shown with solid triangles in Figure 3 . Notice that the data show a decrease with frequency which is greater than 
the (frequency)' 1 decrease characteristic of the data at lower frequencies. The Weidman et. al., (1981) results were 
obtained using only the first few microseconds of the radiation waveform. Since, this portion of the waveform is 
most strongly affected by propagation losses, it is not clear whether the rapid decrease evident in the data is a real 
characteristic of first return strokes or an artifact of the measurements. On the other hand, there is some theoretical 
evidence to suggest that a decrease in the spectrum as (frequency)" 2 should manifest itself at high frequencies (Le 
Vine, 1980). 

A few measurements using the Fourier transform approach have been made of the spectrum of events other 
than return strokes. All of these have been reported recently by Weidman, et. al., (1981). Figure 4 shows the 
spectra obtained for positive (top) and negative (bottom) intracloud events. The curves show the spectra of the in- 
tracloud events superimposed on the spectrum of first return strokes (solid line). On the left in Figure 4 are exam- 
ples of positive and negative intracloud events recorded by the author (Le Vine) at the Goddard Space Flight 
Center, in Greenbelt, Maryland in 1982. Notice that at the high frequencies the spectra for these events tend to 
coincide with the spectrum of return strokes, but at low frequencies the spectrum is smaller than the spectrum of re- 
turn strokes and decreasing. This is the behaviour one would expect of a discharge with the same general physical 
characteristics as a return stroke but of shorter length (e.g. Le Vine, 1980). 


III. SPECTRA FROM DIRECT MEASUREMENT 

In addition to the approach described above, it is also possible to determine the spectrum directly by measur- 
ing the power incident at a particular frequency. The procedure is to use a filter to accept signals only in a narrow 
band of frequencies near the one of interest and then to use a detector to measure the power being radiated in this 
frequency band. The precise relationship between the spectrum at the nominal frequency, v 0 , to which the filter is 
tuned and the output from this system depends on the specific filter and detector employed. An important special 
case occurs when the filter is very narrow (a small percentage of the center frequency, v G ) and the detector is an 
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Figure 3 . Spectrum of first return strokes normalized to 50 km. 
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Figure 4 . Waveforms of intracloud events (left) and the spectrum of such events (right) as reported by W eidman et . al . 
(1981). The solid lines on the right are the spectrum of first return strokes (from Figure 3) shown for reference. 










envelope detector of the type employed in conventional AM radio receivers. In this case, if the input is a single 
impulse, the spectrum S(v D ) is proportional to the peak value, e p , of the output. This relationship is derived in 
Appendix A where it is shown that with an ideal bandpass filter of bandwidth B and system gain G one obtains: 


S(v 0 ) = 


2GB 


( 1 ) 


Data obtained using this technique are illustrated in Figure 5. The amplitude of the spectra have been nor- 
malized to lightning at 50 km using (distance)' 1 as the rule for scaling amplitude and (bandwidth)' 1 as the rule for 
converting from peak field measurements to spectra when necessary (Equation 1). Among the earliest measure- 
ments of this type are those of Schafer and Goodall (1939) who were interested in assessing the potential interfer- 
ence lightning presented for transmission of television signals. Their data at 139 MHz is shown with a solid 
triangle in Figure 5. Extensive measurements over a wide range of frequencies were made by Homer and Bradley 
(1964). These data are shown with the solid line in Figure 5. Additional measurements by Kosarev et. al., (1970) 
and in 1963 by Hallgren and McDonald (adapted from Cianos, Oetzel and Pierce, 1972) are also shown in Figure 
5 . The spectrum which emerges from these measurements is one which peaks near 5 kHz and then decreases as 1/f 
to the limits of the measurements (1 GHz). 


The preceeding are only some of the measurements which have been made with this technique. A large 
number of measurements have been made, but under a great variety of circumstances (e.g., distance, type of re- 
ceiving equipment and bandwidth). It is often difficult to reduce the measurements to common units. Nevertheless, 
several attempts have been made to combine the measurements by normalizing the data to common bandwidth and 
distance (e.g., Homer, 1964,; Oh, 1969; Kimpara, 1965; Cianos, Oetzel and Pierce, 1972). Additional data 
gathered from these reviews are plotted in Figure 6. The relationship (distance)" 1 has been used to convert to light- 
ning 50 km away and the relationship (bandwidth)' 1 has been used to convert peak field measurements to spectra 
(Equation 1). Shown with open circles are data by Takagi and Takeuti (1963) as reported by Kimpara (1965). Data 
by Iwata and Kanada (1967) as reported by Cianos, Oetzel and Pierce (1972) are shown with X’s. Finally, mea- 
surements of several researchers using radar receivers (Atlas, 1959; Hewitt, 1957; Pawsey, 1957) as reported by 
Oh (1969) and Cianos, Oetzel and Pierce (1972) are shown with open triangles. Notice that these additional data 
tend to fall somewhat below the 1/f spectrum suggested in Figure 5 and with a significant spread of the data, espe- 
cially at the higher frequencies. 

A major problem faced in interpretating spectra such as presented in Figures 5 and 6 is that the data do not rep- 
resent radiation from single events (e.g. return strokes) but rather represent some generally unknown collection of 
events in the flash. To illustrate the nature of the problem, the output from several AM radio receivers as seen dur- 
ing a representative cloud-to-ground lightning flash is shown in Figure 7. The records at 3, 30 and 300 MHz are 
data from real lightning recorded by the author in Florida during the Thunderstorm Research Project (TRIP-76; 
Pierce, 1976) using standard AM radio receivers with a common bandwidth of 300 kHz. The records at 30 and 300 
kHz are the author’s impression of what radiation from the flash would look like based on reports in the literature 
(e.g., Homer and Bradley, 1964; Malan, 1958). The slow electric field change for this flash is shown at the bottom 
for reference. Notice that the radiation consists of many discrete impulses. Some of the impulses correlate with 
identifiable portions of the discharge. For example, the initial impulses in this record occur where one expects to 
find the stepped leader, and the return strokes (which occur at the abrupt changes in the slow electrical field change 
record) are associated with large impulses of RF radiation. However, these are only a few of the impulses seen dur- 
ing the flash. The others are probably associated with intracloud portions of the discharge about which we have 
only begun to learn. The discrete nature of the radiation is characteristic of the data recorded by this author from 
both cloud-to-ground and intracloud discharges (e.g., Le Vine, 1976). Another important characteristic of the data 
is that the radiation does not appear to change from a series of discrete impulses at the lower frequencies to a con- 
tinuum at the higher frequencies as reported in some investigations (e.g., Homer and Bradley, 1964; Malan, 
1958). Rather the radiation consists of a sequence of discrete impulses at all frequencies in the range (3 - 300 
MHz) investigated by this author. In addition, the impulses tend to correlate well among the frequencies. That is, a 
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Figure 5 . Spectrum of radiation from lightning obtained using radio frequency (fdter-detector) techniques 
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Figure 6. Spectrum of radiation from lightning obtained using radio frequency techniques. 






30 kHz 



Figure 7. Radiation from a cloud-to-ground lightning flash. The signals at 3-300 MHz are data 
recorded by the author. The examples at 30 and 300 kHz are an estimate based on measurements 
reported in the literature. 


10 




particular impulse tends to be present at all frequencies. The behaviour is consistent with a lightning flash consist- 
ing of a sequence of individual discharges (leader steps, return strokes, K-changes, etc.) each radiating over a very 
broad band of frequencies. Certainly the return stroke is a discharge with this characteristic (Figure 3). 

The problem which occurs with spectral measurements using the filter-detector approach is distinguishing be- 
tween impulses. The impulse response of the measuring system is inversely proportional to the system bandwidth; 
consequently there exists a fundamental contradiction between the requirement for a bandwidth narrow enough to 
be close to the nominal frequency of interest and the ability of the measurement to distinguish between closely 
spaced impulses. For example, with a bandwidth of 250 Hz, the impulse response of the system is on the order of 
0.25 seconds. This would include many impulses in a typical lightning discharge (Figure 7). Very few inves- 
tigators have attempted to distinguish between events with measurements of this type and those who have tried 
have had difficulty (e.g., Takagi and Takeuti, 1963). The spectra which have been reported using this method 
generally must be regarded as an integral (average) over many events in the flash. The exception is at VLF fre- 
quencies where, because of the low frequency, return strokes tend to be the dominant source of radiation. 


IV. DISCUSSION 

Notice the similarity between the spectrum of the first return stroke (Figure 3) and the composite spectrum 
obtained from all the filter-detector measurements (Figures 5 and 6). To facilitate the comparison, the spectra 
have been plotted together in Figures 8 and 9. Figure 8 shows the return stroke spectrum (Figure 3) and the spec- 
trum in Figure 5 together, and Figure 9 shows the return stroke spectrum and all the filter-detector measurements 
(Figure 6) together. Considering all the variables that enter such measurements, the spectra are very similar. This 
is especially so at frequenices below 1 MHz, but even at higher frequencies where the spread in the data is great, 
the two spectra overlap. 

The similarity between the spectra is not surprising because of any fundamental difference between the two 
techniques. In fact, in principle the two techniques for measuring the spectrum of lightning discharges ought to 
yield identical results. In practice, they should be complimentary, the Fourier transform approach having advan- 
tages at low frequencies and the direct (filter-detector) approach having advantages at higher frequencies. How- 
ever, because of the relatively narrow bandwidth employed, the filter-detector technique does not measure radia- 
tion from a single event, but rather is an average over many events in the flash. Furthermore, the data in Figures 
5 and 6 were obtained from measurements of different bandwidth and without any attempt at identifying the por- 
tion of the flash monitored. Consequently, it would seem reasonable to assume that Figures 5 and 6 are more rep- 
resentative of the composite flash rather than of any particular event. What is surprising, then, is that the spec- 
trum of the composite flash and the spectrum of one particular event, the return stroke, are so similar. 

The similarity between the spectrum of one event, the return stroke, on the one hand, and what probably rep- 
resents the spectrum of the composite flash on the other hand, suggests to this author a physical process common 
to all the various individual discharges which make up the lightning flash. There is some additional evidence sup- 
porting this view. For example, Weidman, et. al., (1981) measured spectra of intracloud processes and stepped 
leaders, and these tend to fall on the spectra for return strokes (Figure 4). The difference is a roll- off at the low 
frequencies which is to be expected for identical discharges of shorter length. Also, Le Vine (1976, 1980) has 
been able to predict spectra for the return stroke and composite flash which agree with measurements by assuming 
a transmission line model common to all events. (The events differ in such parameters as channel length, peak 
current and tortuosity.) 

A critical test of this hypothesis occurs at high frequencies where all events should behave similarly. But this 
is where the data is most scattered. Consequently, it would seem important to obtain reliable spectral measure- 
ments of individual events at frequencies above a few MHz. This would resolve the ambiguities in the spectral 
shape apparent in Figures 5 and 6 and would add insight into the physics of the lightning discharge. 
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Figure 8. Composite spectrum of radiation from lightning normalized to 50 km. (Figure 3 and 5 combined). 
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Figure 9. Composite spectrum of radiation from lightning normalized to 50 km. (Figures 3 and 6 combined). 
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APPENDIX A 


RADIO RECEIVER RESPONSE. 


It is the purpose of this appendix to derive an expression for the impulse response of an RF receiving system 
typical of the sort used to measure the spectrum of radiation from lightning. The system consists of an antenna, 
usually vertically polarized, connected to a standard AM radio receiver and followed with a post detection filter 
(Figure Al). Assuming that the antennas are vertically polarized and isotropic over the bandwidth of the measure- 
ments, they can be assumed to deliver a voltage proportional to the incident (vertical) electric field. Except for a 
phase which is ultimately lost in the detector, this proportionally constant is V uG/A e /k where k = wVlIi = 2tt/ 
A., G is the gain of the antenna in the plane parallel to the ground, and A e is the effective receiving area of the an- 
tenna (Kraus, 1966). The receiver is a device for detecting and amplifying the envelope of an amplitude mod- 
ulated sinusoid (carrier) at a particular frequency. This is normally done by translating the input signal to an inter- 
mediate frequency where the actual processing is done. However, the frequency translation is done for engineer- 
ing purposes to make the detection more efficient, and it is not necessary to do this in order to model the receiver 
output. The ideal device is a perfect envelope detector in series with a filter which represents the equivalent 
bandwidth and gain of the system. 


For systems whose bandwidth 8 is small compared to the nominal frequency, v c , of the measurement, these 
operations can be written explicitly in terms of the Fourier transform of the incident radiation. To do so it is con- 
venient to write the vertical component of incident radiation, E(t), in the form: 


E(t) = 


r 

J o 


Re I 2E (v)e' j2lTVt dv 


(Al) 


where Re means “real part of’ and where E(v) is the Fourier transform of E (t). The integral is called the complex 
analytic representation of E (t) (e.g.. Bom and Wolf, 1959). Using this notation, the signal V(t) out of the antenna 
and filter and into the envelope detector is: 

J *O0 

2a(v)H(v)E (v)e' j2m ' t dv (A2) 

0 


where H(v) is the Fourier transform of the filter h(t) and a(v) = (A/k) V^G/Ag is the combined effect of the. an- 
tenna and an amplifier with gain, A. Although in practice the post detection filter is applied to the video output of 
the receiver, mathematically the effect of the post detection filter can also be included in H(v), This will be done 
here and it will be assumed that the equivalent filter H(v) has a passband centered about frequency v c which is 
very narrow compared to v 0 . That is: 

{ H(v - v c ) |v-vj<8 

(A3) 

O |v-v 0 |>8 


and 8/v c << 1 . Using this assumption and making the change of variables £ = v-v 0 , Equation A2 becomes 


V(t) = Re 



2a(v c + £) E (v D + £) H(£)e' j2ll€t d£ 


} 


= j e (v 0 ,t) j cos(2irv 0 t + <J>) 


(A4) 


A-l 



where 


e (v 0 ,t) 



2a(v 0 + 0E ( Vo + 0H(0 e- j2lt€t dC 


(A5) 


and where the assumption that H(v) is nonzero only in a narrow band about v Q (Equation A3) has been used to for- 
mally extend the lower limit of the integration in Equation A4 to infinity. From Equation A4 it is clear that V(t) has 
the form of an amplitude modulated sinusoid, cos(2irv 0 t+4>), at frequency v 0 . The output of the detector is the en- 
velope, |e (v o; t)|, of this carrier. Thus, the output, e 0 (t) can be written: 


e c (t) = 



2a(v o +0E(v o + 0H(0e- j2 ^‘ 



(A6) 


Now, suppose that over the narrow band of frequencies passed by the filter a(v D + 0 = a(v c ) and E(v 0 + 0 = 
E(v 0 ). Then, 


e 0 2 (t) = |E(v 0 )| 2 h 2 (t) 

where |E(v e )| is the magnitude of the spectrum of the electric field at frequency v c and where 

h(t)= I °° 2 a(v o )H(0 e" j2lT ^ dij 

J -00 


(A7) 


(A8) 


is the impulse response of the system. (Note: e 2 = e e*) Integrating both sides of Equation A7 and using Parseval’s 
theorem for Fourier transforms, one obtains: 



= 4|E(v o )| 2 J'"|a(v o )H(0| 2 d4 


Finally, solving for the spectrum, one obtains: 


where 



(A9) 


(A10) 


(All) 


The integral in Equation A1 1 is commonly called the bandwidth (power bandwidth) of the system and |a(v G ) H(0)| 
is the gain of the system. Equation A 10 states that the spectrum of the input signal can be obtained by integrating 
the output and dividing by twice the gain times the square root of the bandwidth of the system. This result only 


A-2 



applies to individual input signals, E(t), whose bandwidth is much larger than the bandwidth of the system (i.e. 
appear as impulses on the time scale of the impulse response of the system). This restriction is a consequence 
of factoring E(v 0 ) out of Equation A6. If the input were a sequence of such pulses rather than an individual event 
E(v) in Equation A6 would be a sum of the form £ E s (v) e’ 2lTvt ' where E;(v) is the spectrum of the individual 
pulses and t ; is the time between pulses. In this case the spacing between pulses can affect the spectral estimate 
(Le Vine, 1977; Dennis and Pierce, 1964). 


An important special case occurs when the system can be modelled as an ideal bandpass filter: 

Itl < B/2 
l£l > B/2 

Then, from Equation A6 one obtains: 

* B/2 


H(0 = 


r h g 

l o 


(A 12) 


e o(0 = 2|a(v 0 ) E(v c ) H 0 J exp(-j2irvt) dv| 
J -B/2 

= 2|a(v 0 )E(v 0 )|H 0 Bsinc(TTBt) 

= e p sine (irBt) 


(A13) 


where e p is the peak value of the signal out of the receiver. Now, squaring and integrating over all time, one obtains: 


/ eo(t) dt = e p J sine 2 (irBt) dt 

00 J _00 


= e p /B 


(A 14) 


Finally, putting this result into Equation A10, one obtains: 


|E(v 0 )| = 


2GB 


(A 15) 


where G = |a(v c ) H(0) | is the gain of the system. Equation A15 applies only if the time between pulses is long 
compared to the response time of the system. 

Another important special case occurs when the input is a random process. If the process consists of a sequ- 
ence of identical pulses with random amplitude and/or arrival time, then the analysis proceeds as above with mod- 
ifications as indicated in the text and described in detail by Le Vine (1977). However, if the random process is 
noise-like (i.e. a continuous, fluctuating signal), then the analysis must be modified. In this case, the appropriate 
definition of the spectrum, S(v), is the Fourier transform of the autocorrelation function, R(t) = <E(t)E* (t + 
t)> of the input signal. Assuming that the detector is an id eal squ are-law detector and that the input, E(t), is a 
stationary, ergodic random process, and defining |E(v)| = V|S(v)|, one obtains 


|E(v 0 )| = 



(A 16) 


\I/Kp rp A ic oc in All 

--Ilvlv Li lo uo Uviiuvvi in i^v|uauuil ill A 


the case of an ideal bandpass filter (Equation A 12), one has A = ja(v Q ) 
H(0) | 2 B. In this case. Equation A16 shows that jE(v 0 )) is proportional to 1/VB and can be obtained from the RMS 
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value of the output signal by dividing by the square root of the bandwidth and twice the system gain. This result 
has been employed to compute spectra of lightning in some cases (e.g. Oh, 1969); however, as indicated in the 
text, lightning is intrinsically impulsive in nature and as a result this formula must be used with caution. 



Figure A1 . Example of a radio receiver used to measure spectrum of radiation from lightning. 
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